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ABSTRACT 

We have analysed the mass and velocity distributions of two samples of relaxed 
elliptical-like-objects (ELOs) identified, at z — 0, in a set of self-consistent hydro- 
dynamical simulations operating in the context of a concordance cosmological model. 
ELOs have been identified as those virtual galaxies having a prominent, dynamically 
relaxed stellar spheroidal component, with no extended discs and very low gas content. 
Our analysis shows that they are embedded in extended, massive dark matter haloes, 
and they also have an extended corona of hot diffuse gas. Dark matter haloes have ex- 
perienced adiabatic contraction along their assembly process. The relative ELO dark- 
to bright-mass content and space distributions show broken homology, and they are 
consistent with observational results on the dark matter fraction at the central regions, 
as well as on the gradients of the mass-to-light ratio profiles for boxy ellipticals, as 
a function of their stellar masses. These results indicate that massive ellipticals miss 
stars (i.e., baryons) at their central regions, as compared with less massive ones. Our 
simulations indicate that these missing baryons could be found beyond the virial radii 
as a hot, diffuse plasma. This mass homology breaking could have important implica- 
tions to explain the physical origin of the Fundamental Plane relation. The projected 
stellar mass profiles of our virtual ellipticals can be well fit by the Sersic formula, with 
shape parameters n that agree, once a stellar mass-to-light ratio independent of posi- 
tion is assumed, with those obtained from surface brightness profiles of ellipticals. The 
agreement includes the empirical correlations of n with size, luminosity and velocity 
dispersion. The total mass density profiles show a power-law behaviour over a large 
p/^Vir interval, consistent with data on massive lens ellipticals at shorter radii. The 
velocity dispersion profiles show kinematical segregation, with no systematic mass de- 
pendence (i.e., no dynamical homology breaking) and a positive anisotropy (i.e., radial 
orbits), roughly independent of the radial distance outside the central regions. The 
LOS velocity dispersion profiles are declining. These results give, for the first time 
from cosmological simulations, a rather detailed insight into the intrinsic mass and 
velocity distributions of the dark, stellar and gaseous components of virtual ellipti- 
cals. The consistency with observations strongly suggests that they could also describe 
important intrinsic characteristics of real ellipticals, as well as some of their properties 
recently inferred from observational data (for example, downsizing, the appearance of 
blue cores, the increase of the stellar mass contributed by the elliptical population as 
z decreases). 

Key words: galaxies: elliptical and lenticular, cD - galaxies: haloes - galaxies: kine- 
matics and dynamics - galaxies: structure - dark matter - hydrodynamics 
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Among all galaxy families, elliptical galaxies (EGs) are the 
simplest ones and those that show the most precise empiri- 
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cal regularities, such as: the very similar shapes their surface 
brightness profiles show (i.e., the Sersic law, Sersic f968), 
the power-law form of their three-dimensional mass density 
profiles emerging from some strong lensing analyses (Koop- 
mans et al. 2006), and the relations among some of their 
observable parameters. The importance of these regularities 
lies in that they very likely encode a lot of relevant infor- 
mations on the physical processes underlying the histories 
of the mass assembly, the rates of dissipation and the rates 
of star formation (SF) of EGs. 

Despite of their interest, very few is known, both from 
the theoretical or observational points of view, about the 
mass or velocity distributions of the different ellipticals mass 
components (stars, hot and cold gas and dark matter). There 
has been, nevertheless, an important recent progress on the 
photometric characterisation of EGs, and, in fact, authors 
now agree that the Sersic law adequately describes the op- 
tical surface brightness profiles of most of them (Caon, Ca- 
paccioli, & D'Onofrio 1993; Trujillo, Graham, & Caon 2001; 
Bcrtin, Ciotti, & Del Principe 2002). The Sersic law can be 
written 

jlight (jR) = /l^t e;cp[ _ 6n(fl/Ji ligh t) l/„ ]; (1) 

where I hsht (R) is the surface brightness at projected 
distance R from the ellipticals centre, R]l sht is the effective 
half-light radius, encompassing half the total galaxy lumi- 
nosity, b n ~ 2n — 1/3 + 0. 009876/n, and n is the Sersic shape 
parameter. Putting n — 4 the largely used de Vaucouleurs 
R 1 / 4 law (de Vaucouleurs 1948) is recovered. 

It is generally assumed that galaxies of any type are 
embedded in massive haloes of dark matter. However, from 
the observational point of view, the importance and the dis- 
tribution of dark matter in EGs is still a matter of a liv- 
ing debate. Data on stellar kinematics from integrated-light 
spectra are very scarce beyond 2i?e ght , making it difficult 
even to establish the presence of a dark matter halo (Kron- 
awitter et al. 2000; Magorrian & Ballantyne 2001) through 
this method. Otherwise, the lack of mass tracers at larger 
distances that can be interpreted without any ambiguity, has 
historically hampered the proper mapping of the mass dis- 
tribution at the outer regions of EGs. The situation is chang- 
ing and a dramatic improvement is expected in the near fu- 
ture. In fact, several ongoing projects have already produced 
high quality data on samples of ellipticals through differ- 
ent methods, for example: stellar kinematics from integral- 
field spectroscopic measurements (SAURON; de Zeeuw et 
al. 2002; Cappellari et al. 2006); strong gravitational lens- 
ing (CLASS; Myers et al. 1995; LSD; Koopmans & Treu 
2003, Treu & Koopmans 2004; SLACS; Koopmans et al. 
2006); stellar kinematics from planetary nebulae (PNs; Dou- 
glas et al. 2002), or globular cluster (Bergond et al. 2006) 
observations; and X-rays (O'Sullivan & Ponman 2004a,b). 
In particular, the combination of high-quality stellar spec- 
troscopy and strong lensing analyses breaks the so-called 
mass-anisotropy degeneracy, giving strong indications that 
constant mass-to-light ratios can be ruled out at > 99% 
CL, consistent with the presence of massive and extended 
dark matter haloes around, at least, the massive lens ellip- 
ticals analysed so far (Treu & Koopmans 2004; Koopmans 
et al. 2006). Moreover, these authors have also found that 
the dark matter and the baryons mass density profiles com- 
bine in such a way that the total mass density profiles can 



be fit by power-law expressions within their Einstein radii, 
whose slopes are consistent with isothermality. Similar con- 
clusions on the important amounts of dark matter inside the 
virial radii of ellipticals have been reached from weak lens- 
ing of L, galaxies (Guzik & Seljak 2002; Hoekstra, Yee & 
Gladders 2004), dynamical satellite studies (van den Bosch 
et al. 2004) and X-ray analyses (Humphrey et al. 2006). 
Other observational results or some of their interpretations, 
however, could suggest that the amounts of dark matter in 
the haloes of some ellipticals are not that important. For 
example, Napolitano et al. (2005) have analysed the mass- 
to-light gradients of a sample of elliptical + SO galaxies, 
and found that these are positive and important in massive, 
boxy EGs, but no very important for faint, discy EGs. This 
has been confirmed by Ferreras, Saha & Williams 2005 using 
lensing analyses. This result is similar to what Romanowsky 
et al. (2003; see also Romanowsky 2006) have propounded 
from the study of random velocities at the outskirts of EGs 
through PN, found to be low, and first interpreted by these 
authors as proving a dearth of dark matter in EGs, while 
Dekel et al. (2005) explain these large-radii low velocity dis- 
persions as an effect of anisotropy and triaxiality of the halo 
stellar populations of these galaxies. 

Assuming that ellipticals are embedded in massive 
haloes of dark matter, a second important concern is the 
possibility that their profiles have near-universal shapes. 
Here most inputs come from numerical simulations because 
observational inputs are scarce. When no dissipative pro- 
cesses are taken into account, spherically averaged dark mat- 
ter density profiles of relaxed haloes produced in N-body 
simulations have been found to be well fitted by analyti- 
cal expressions such that, once rescaled, give essentially a 
unique mass density profile, determined by two parameters. 
These two parameters are usually taken to be the total mass, 
M v i r , and the concentration, c or the energy content, E. 
These two parameters are, on their turn, correlated (i.e., the 
mass-concentration relation, see, for example, Bullock et al. 
2001; Wechsler et al. 2002; Manrique et al. 2003). When hy- 
drodynamical forces and cooling processes enter the assem- 
bly of these haloes and the baryonic objects they host, the 
dark matter profiles could be modified in the regions where 
baryons are dynamically dominant, due the so-called adia- 
batic contraction (see, for example, Blumenthal et al. 1986; 
Dalcanton, Spergel, & Summers 1997; Tissera & Dommguez- 
Tenreiro 1998; Gnedin et al. 2004; Gustafsson, Fairbain & 
Sommer-Larsen 2006). So, the shapes of dark matter haloes 
in ellipticals could deviate from the near-universal behaviour 
of dark matter haloes produced in purely N-body simula- 
tions. 

Another important issue concerns the three dimensional 
cold baryon mass (i.e., stellar mass) distribution, and, more 
particularly, its distribution relative to the dark matter 
haloes: are ellipticals homologous systems or is the homology 
broken in their stellar mass distribution or in their relative 
dark- versus bright-mass distribution? 

Concerning the three dimensional velocity distributions 
of the different elliptical components, very few is known 
either. In particular, the anisotropy of the stellar three- 
dimensional velocity dispersion tensor is hard to be observa- 
tionally characterised. This is an important issue, however, 
not only because anisotropy plays an important role in the 
analyses of the elliptical dark matter content at several effec- 
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tive radii, but also because it could keep fossil informations 
about the physical processes involved in mass assembly and 
stellar formation in EGs. The relative behaviour of the three- 
dimensional velocity dispersion tensors for the stellar and 
the dark mass components (i.e., the so-called kinematical 
segregation) is still more uncertain. There is not an unam- 
biguous observational input about its presence in ellipticals, 
or about its possible systematic dependence with the ellip- 
tical mass scale. This is an important point because of its 
possible connection with the physical origin of the so-called 
Fundamental Plane (FP) relation, as different authors have 
suggested (Graham & Colless 1997; Busarello et al. 1997; 
Pahre, de Carvalho & Djorgovski 1998). 

Understanding the FP relation (Djorgovski & Davis 
1987; Dressier et al. 1987; Faber et al. 1987; Bernardi et 
al. 2003 and references therein) is in fact a milestone in 
the physics of elliptical formation. The FP is defined by the 
three observational parameters characterising the mass and 
velocity distributions of the stellar component: luminosity L, 
radius at half projected light i?e ght i an d central line-of-sight 
velocity dispersion <7i OSi o- The FP is tilted relative to the 
virial plane. Different authors interpret this tilt as caused 
by different misassumptions relative to the constancy of the 
dynamical mass-to-light ratios, M v i r /L, or the mass struc- 
ture coefficients, cJ,J r = 2 GMvi 1 r ight , with the mass scale (see 

discussion in Onorbe et al. 2005, 2006, and references cited 
therein). Recently, Onorbe et al. (2005, 2006) have found 
that the samples of elliptical-like objects identified at z — 
in their fully-consistent cosmological hydrodynamical simu- 
lations show systematic trends with the mass scale in both, 
the relative global dark-to-bright mass content, as well as in 
the relative sizes of the stellar and the dark mass compo- 
nents, giving rise to dynamical FPs. 

We see that the mass or velocity distributions of the 
different elliptical mass components encode a lot of infor- 
mations about the physical origin of the FP, and, conse- 
quently, on the physics of elliptical formation. We see also 
that, unfortunately, observational methods, by themselves, 
suffer from some drawbacks to deepen into these issues. A 
major drawback is that the informations on the intrinsi- 
cal mass distribution are not directly available: we see the 
projections of luminosity (not three-dimensional mass) dis- 
tributions. Another major drawback is that the intrinsic 3D 
velocity distribution of galaxies is severely limited by pro- 
jection, only the line-of-sight velocity distributions can be 
inferred from galaxy spectra. And, so, the interpretation of 
observational data is not always straightforward. To com- 
plement the informations provided by data and circumvent 
these drawbacks, analytical modelling is largely used in lit- 
erature (Kronawitter et al. 2000; Gerhard et al. 2001; Ro- 
manowsky & Kochanek 2001; Borriello, Salucci & Danese 
2003; Padmanabhan et al. 2004; Mamon & Lockas 2005a, 
2005b) . They give very interesting insights into mass and ve- 
locity distributions, as well as the physical processes causing 
them, but are somewhat limited by symmetry considerations 
and other necessary simplifying hypotheses. These difficul- 
ties and limitations could be circumvented should we have 
at our disposal complete informations on the phase-space 
of the galaxy constituents. This is not possible through ob- 
servations, but can be attained, at least in a virtual sense, 
through numerical simulations. 



The first authors who studied the formation and proper- 
ties of EGs by means of numerical methods used purely grav- 
itatory pre-prepared simulations. The origin of the FP was 
addressed, among others, by Capelato, de Carvalho & Carl- 
berg (1995); Gonzalez- Garcia & van Albada (2003); Dantas 
et al. (2003); Nipoti, Londrillo & Ciotti (2003, 2006) and 
Boylan-Kolchin, Ma & Quataert (2005). Some of these au- 
thors (Dantas et al. 2003, Nipoti et al. 2003) conclude that 
dissipation must be a basic ingredient in elliptical forma- 
tion. Bckki (1998) first considered the role of dissipation 
in elliptical formation through pre-prepared simulations. He 
adopts the merger hypothesis (i.e., ellipticals form by the 
mergers of two equal-mass gas-rich spirals) and he focuses 
on the role of the timescale for SF in determining the struc- 
tural and kinematical properties of the merger remnants. 
He concludes that the slope of the FP reflects the differ- 
ence in the amount of dissipation the merger end products 
have experienced according with their luminosity (or mass) . 
Recently, Robertson et al. (2006) have confirmed this con- 
clusion on the role of dissipative dynamics to shape the FP, 
again through pre-prepared mergers of disc galaxies. 

Apart from the FP relation origin, other aspects of the 
formation and evolution of EGs have been analysed through 
pre-prepared simulations. For example, the different isopho- 
tal shapes of ellipticals and their kinematical support (either 
intrinsic rotation or anisotropic dispersion) have been ad- 
dressed by A. Burkert and coworkers, who analysed in detail 
binary mergers of virtual galaxies with different morpholo- 
gies and initial conditions (Naab & Burkert 2003; Naab & 
Trujillo 2006; Naab, Khochfar & Burkert 2006). Otherwise, 
Gonzalez-Garcfa, Balcells & Olshevsky (2006) analysed the 
velocity distribution of dissipationless binary merger rem- 
nants involving galaxies with different morphologies. 

We see that pre-prepared simulations of merger events 
provide a very useful tool to work out the mass and veloc- 
ity distributions of EGs. They allow also to find out their 
links with the processes involved in galaxy assembly, but 
they are somewhat limited, for example by the fact that the 
probability of a particular initial setup at a given z is not 
known a priori, and that mergers involving more than two 
objects also occur and are frequent at high zs, so that some 
complementary informations must be provided, for example 
through semi-analytical models (Khochfar & Burkert 2005; 
Naab, Khochfar & Burkert 2006). 

To overcome these limitations, a convenient method is 
to study the processes involved in galaxy formation in a 
cosmological context through self- consistent gravodynamical 
simulations. The method works as follows: initial conditions 
are set at high z, in an homogeneously sampled box, as 
a Montecarlo realization of the field of primordial fluctu- 
ations to a given cosmological model; then the evolution of 
these fluctuations is numerically followed up to z = by 
means of a computer code that solves, in a periodic box, 
the gravitational and hydrodynamical evolution equations. 
This method allows us to follow the evolution of the dynam- 
ical and thermohydrodynamical properties of matter in the 
universe; individual galaxy-like objects naturally appear as a 
consequence of this evolution. No prescriptions are needed as 
far as their mass assembly processes are concerned, only SF 
processes need further modelling. The important point here 
is that self-consistent simulations directly provide with com- 
plete 6-dimensional phase-space informations on each con- 
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stituent particle sampling a given galaxy-like object formed 
in the simulation, that is, they give directly the mass and 
velocity distributions of dark matter, gas and stars of each 
object. 

Kobayashi (2005) has simulated the chemodynamical 
evolution of 74 fields with different cosmological cold dark 
matter initial spectra set in slowly rotating spheres, each of 
them with a 1.5 Mpc comoving radius and vacuum bound- 
aries. So, these simulations are not yet fully self-consistent. 
She succeeded in reproducing the observed global scaling re- 
lations shown by EGs, and, in particular, the FP relation, 
and the surface-brightness profiles, as well as the colour- 
magnitude and the mass-metallicity relations. She also anal- 
yses the role of major merger events and the timescales for 
SF in shaping the mass and sizes of remnants. 

Concerning self-consistent hydrodynamical simulations, 
Sommer-Larsen, Gotz & Portinari (2002) present first re- 
sults on early-type galaxy formation in a cosmological con- 
text. Meza et al. (2003) present results of the dissipative for- 
mation of a compact elliptical galaxy in the ACDM scenario. 
Kawata & Gibson (2003, 2005) and Gibson et al. (2006) have 
studied the X-ray and optical properties of virtual ellipticals 
formed in different simulations run with their chemodynam- 
ical Tree/SPH code. Romeo, Portinari & Sommer-Larsen 
(2005) analyse the galaxy stellar populations formed in their 
simulations of galaxy clusters. Naab et al. (2005) got, from 
cosmological initial conditions, a spheroidal system whose 
photometric and kinematical properties agree with observa- 
tions of ellipticals, in a scenario not including feedback from 
supernovae or AGN and not requiring recent major mergers. 
Interesting results on elliptical formation have also been ob- 
tained by de Lucia et al. (2006), from a semi-analytic model 
of galaxy formation grafted to the Millennium Simulation. 

However, detailed analyses of the mass and velocity dis- 
tributions of samples of virtual ellipticals formed in fully 
self-consistent hydrodynamical simulations, and, in particu- 
lar, of the amount and distribution of dark matter relative 
to the bright matter distribution, as well as of the kinemat- 
ics of the dark and bright components, and their successful 
comparison with observational data, were still missing. This 
is the work we present in this paper. To be specific, we have 
analysed the samples of virtual ellipticals formed in ten self- 
consistent hydrodynamical simulations, run in the frame- 
work of a flat ACDM cosmological model characterised by 
cosmological parameters consistent with their last determi- 
nations (Spergel et al. 2006). Galaxy- like objects of differ- 
ent morphologies appear in these simulations at z = 0: disc- 
like objects, SO-like objects, elliptical-like objects (hereafter, 
ELOs) and peculiars. 

Our previous work has shown that these ELO samples 
have properties that agree with observational data, so that 
they have counterparts in the real world. In fact, an analysis 
of the structural and dynamical ELO parameters that can 
be constrained from observations (i.e., stellar masses, pro- 
jected half-mass radii, central line-of-sight velocity disper- 
sions), has shown that they are consistent with those mea- 
sured in the SDSS DR1 elliptical sample (Saiz, Dommguez- 
Tenreiro & Serna 2004), including the FP relation (Onorbe 
et al. 2005), and their lack of evolution at low and interme- 
diate redshifts (Dommguez-Tenreiro et al. 2006, hereafter 
DTal06). Also, their stellar populations have age distribu- 
tions showing similar trends as those inferred from obser- 



vations, i.e., most stars have formed at high z on short 
timescales, and, moreover more massive objects have older 
means and narrower spreads in their stellar age distribu- 
tions than less massive ones * (Dommguez-Tenreiro, Saiz & 
Serna, 2004, hereafter DSS04). 

The paper is organised as follows: in §2 we briefly de- 
scribe the simulations and the SF modelling. In §3, the ELO 
sample and the generic structure of ELOs are described. A 
brief account on ELO formation is given in §4. The three 
dimensional orbital structure of the baryonic component is 
briefly described in §5. §6 is devoted to report on the dark 
matter and baryonic mass density profiles of the ELOs pro- 
duced in the simulations. In §7 we report on the kinematics 
of the different ELO components. Finally, in §8 we sum- 
marise our results and discuss them in the context of theoret- 
ical results on halo structure and dissipation of the gaseous 
component and their connections with observational data. 



2 THE SIMULATIONS 

We have analysed ELOs identified in ten self-consistent cos- 
mological simulations run in the framework of the same 
global flat ACDM cosmological model, with h = 0.65, 
f2 m = 0.35, Qb = 0.06. The normalisation parameter has 
been taken slightly high, erg = 1.18, as compared with the 
average fluctuations of 2dFGRS or SDSS galaxies (Spergel et 
al. 2006) to mimic an active region of the Universe (Evrard, 
Silk & Szalay 1990). 

We have used a lagrangian code (DEVA; Serna, 
Dommguez-Tenreiro & Saiz 2003), particularly designed to 
study galaxy assembly in a cosmological context. Gravity is 
computed through an AP3M-like method, based on Couch- 
man (1991). Hydrodynamics is computed through a SPH 
technique where special attention has been paid that the im- 
plementation of conservation laws (energy, entropy and an- 
gular momentum) is as accurate as possible (see Serna et al. 
2003 for details, in particular for a discussion on the implica- 
tions of violating some conservation laws). Entropy conser- 
vation is assured by taking into consideration the space vari- 
ation of the smoothing length (i.e., the so-called X7h terms). 
Time steps are individual for particles (to save CPU time, 
allowing a good time resolution), as well as masses. Time in- 
tegration uses a PEC scheme. In any run, an homogeneously 
sampled periodic box of 10 Mpc side has been employed and 
64 3 dark matter and 64 3 baryonic particles, with a mass of 
m dark = 1 2g x 1Q 8 and m b = 2 67 x 10 7 Mqj respe ctively, 

have been used. The gravitational softening used was e = 2.3 
kpc. The cooling function is that from Tucker (1975) and 
Bond et al. (1984) for an optically thin primordial mixture 
of H and He (X — 0.76, Y — 0.24) in collisional equilibrium 
and in absence of any significant background radiation field 
with a primordial gas composition. Each of the ten simula- 
tions started at a redshift z; n = 20. 

SF processes have been included through a simple phe- 
nomenological parametrization, as that first used by Katz 
(1992, see also Tissera, Lambas & Abadi 1997 and Serna et 
al. 2003 for details) that transforms cold locally-collapsing 

* this is equivalent to the downsizing concept introduced by 
Cowie et al. (1996) to mean that SF is stronger at low redshifts 
in smaller galaxies than in larger ones 
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gas at the scales the code resolves, denser than a thresh- 
old density, pthrcs, into stars at a rate dp B ta,r/dt — c*p gas /i 9 , 
where t g is a characteristic time-scale chosen to be equal 
to the maximum of the local gas-dynamical time, td yn — 
(47rGp g as) -1 ^ 2 , and the local cooling time; c* is the aver- 
age SF efficiency at resolution e scales, i.e., the empirical 
Kennicutt- Schmidt law (Kennicutt 1998). It is worth not- 
ing that, in the context of the new sequential multi-scale SF 
scenarios (Vazquez- Semadeni 2004a and 2004b; Ballesteros- 
Paredes et al. 2006 and references therein), it has been ar- 
gued that this law, and particularly so the low c* values in- 
ferred from observations, can be explained as a result of SF 
processes acting on dense molecular cloud core scales when 
conveniently averaged on disc scales (Elmegreen 2002, Sar- 
son et al. 2004, see below). Supernova feedback effects or en- 
ergy inputs other than gravitational have not been explicitly 
included in these simulations. We note that the role of dis- 
crete stellar energy sources at the scales resolved in this work 
is not yet clear, as some authors argue that stellar energy 
releases drive gas density structuration locally at sub-kpc 
scales (Elmegreen 2002). In fact, recent MHD simulations of 
self-regulating Type II supernova heating in the interstellar 
medium at scales < 250 pc (Sarson et al. 2004), indicate 
that this process produces a Kennicutt-Schmidt-like law on 
average. If this were the case, the Kennicutt-Schmidt law 
implemented in our code would already implicitly account 
for the effects stellar self-regulation has on the scales our 
code resolves, and our ignorance on sub-kpc scale processes 
would be contained in the particular values of pthrcs and c* . 

Five out of the ten simulations (the SF-A type simula- 
tions) share the SF parameters (pthrcs = 6 x 10~ 25 gr cm~ 3 , 
c» = 0.3) and differ in the seed used to build up the initial 
conditions. To test the role of SF parameterisation, the same 
initial conditions have been run with different SF parame- 
ters (pthrcs = 1.8 x 10~ 24 gr cm~ 3 , c = 0.1) making SF 
more difficult, contributing another set of five simulations 
(hereafter, the SF-B type simulations). 



3 SIMULATED ELLIPTICAL-LIKE OBJECTS 
AT Z=0: THE ELO SAMPLE 

ELOs have been identified as those galaxy-like-objects hav- 
ing a prominent, relaxed spheroidal component, made out of 
stars, with no extended discs and very low cold gas content. 
Moreover, their dark-matter haloes must also be relaxed to 
allow their meaningful analysis, excluding systems that have 
just merged or that are about to merge. It turns out that, 
at z = 0, 26 (17) objects out of the more massive formed in 
SF-A (SF-B) type simulations fulfil these conditions. ELOs 
form two samples (the SF-A and SF-B ELO samples) par- 
tially analysed in Saiz et al. (2004), in DSS04 and in DTal06. 
In Onorbe et al. (2005, 2006) it is shown that both samples 
satisfy dynamical FP relations. 

ELOs are embedded in dark matter haloes whose virial 
radii t are in the range 527 > r v ; r > 191 kpc. ELO stellar 



' The virial radii, r v ; r , have been calculated using the Bryan & 
Norman (1998) fitting formula, that yields, at z = 0, a value of 
A ~ 100 for the mean density within r v ; r in units of the critical 
density 



components have ellipsoidal shapes and have a lower limit in 
their stellar mass content of 3.8 x 10 10 Mq (see Kauffmann 
et al. 2003 for a similar result in SDSS galaxies). Inner discs, 
when present, are made out of cold gas. ELOs have also 
hot diffuse gas forming an extended halo (Saiz, Dommguez- 
Tenreiro & Serna 2003). The number of dark and baryonic 
particles within r v i r are in the ranges (5.3 x 10 4 ,2.4 x 10 3 ) 
and (3.1 x 10 4 , 2.0 x 10 3 ), respectively, giving a lower limit 
in the virial masses of ELOs of M vir > 3.7 x 10 11 M . Some 
ELOs show a clear net rotation, resulting in an average value 
of their spin parameter of A = 0.033. ELO mass function is 
consistent with that of a small group, that is a dense, en- 
vironment (Cuesta-Bolao & Serna, private communication). 
ELOs in the SF-B sample tend to be of later type than their 
corresponding SF-A counterparts because forming stars be- 
comes more difficult; this is why many of the SF-B sample 
counterparts of the less massive ELOs in SF-A sample do 
not satisfy the selection criteria, and the SF-B sample has a 
lower number of ELOs that the SF-A sample. 



4 A BRIEF ACCOUNT ON ELO FORMATION 

The simulations unveil the physical patterns of ELO mass 
assembly, energy dissipation and SF rate histories (see Sierra 
Glez. de Buitrago et al. 2003; DSS04, DTal06). Our simu- 
lations indicate that ELOs are assembled out of the mass 
elements that at high z are enclosed by those overdense re- 
gions R whose local coalescence length L c (t,R) (Vergassola 
et al. 1994) grows much faster than average, and whose mass 
scale (total mass enclosed by R, Mr) is of the order of an 
elliptical galaxy virial mass. Analytical models, as well as N- 
body simulations indicate that two different phases operate 
along halo mass assembly: first, a violent fast one, where the 
mass aggregation rates are high, and then, a slower one, with 
lower mass aggregation rates (Wechsler et al. 2002; Zhao et 
al. 2003; Salvador-Sole, Manrique, & Solanes 2005). Our hy- 
drodynamical simulations indicate that the fast phase occurs 
through a multiclump collapse (see Thomas, Greggio & Ben- 
der 1999) ensuing turnaround of the overdense regions, and 
it is characterised by the fast head-on (that is, with very 
low relative orbital angular momentum) fusions experienced 
by the nodes of the cellular structure these regions enclose, 
resulting in strong shocks and high cooling rates of their 
gaseous component, and, at the same time, in strong and 
very fast star SF bursts that transform most of the avail- 
able cold gas in R into stars. For the massive ELOs in this 
work, this happens between z ~ 6 and ~ 2.5 and mainly 
corresponds to a cold mode of gas aggregation, as in Keres 
et al. (2005) . Consequently, most of the dissipation involved 
in the mass assembly of a given ELO occurs in this violent 
early phase at high z; moreover, its rate history t is reflected 
by the SF rate history, as illustrated in figure 1 of DTal06. 
The Fundamental Plane relation shown by EGs appears in 
this fast violent phase as a consequence of dissipation and 
homology breaking in the mass distribution (see Onorbe et 
al. 2005, 2006 and DTal06). 

The slow phase comes after the fast phase. In this phase, 

t That is, the amount of cooling per time unit experienced by 
those gas particles that at z = form the ELO stellar component 
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the halo mass aggregation rate is low and the M v j r increment 
results from major mergers, minor mergers or continuous 
mass accretion. Our cosmological simulations show that the 
fusion rates are generally low and that these mergers gen- 
erally imply only a modest amount of energy dissipation or 
SF. In fact, a strong SF burst and dissipation follow a ma- 
jor merger only if enough gas is still available after the early 
violent phase. This is unlikely in any case, and it becomes 
more and more unlikely as M v i r increases (see DSS04). A 
consequence of the lack of dissipation is that the Funda- 
mental Plane is roughly preserved along the slow phase (see 
DTal06). We have to point out that mergers play an impor- 
tant role in this slow phase as far as mass assembly is con- 
cerned: an ~ 50% of ELOs in the sample have experienced a 
major merger event at 2 < z < 0, that result in the increase 
of the ELO mass content, size and stellar mean square ve- 
locity. Some of these mergers are multiple, and in some few 
cases, either binary or multiple mergers involve disc galax- 
ies, but none of the ELOs in our sample has been shaped by 
a merger of two or more adult disc virtual galaxies, maybe 
because our simulated box mimics a dense environment. 

So, our simulations suggest that most of the stars of 
to-day ellipticals, or at least of those in dense environments, 
could have formed at high redshifts, while they are assem- 
bled later on [see de Lucia et al. (2005), for similar results 
from a semi-analytic model of galaxy formation grafted to 
the Millennium Simulation]. This formation scenario shares 
some aspects of both, the hierarchical merging and the 
monolithic collapse scenarios, but it has also significant dif- 
ferences, mainly that most stars belonging to EGs form out 
of cold gas that had never been shock heated at the halo 
virial temperature and then formed a disc, as the conven- 
tional recipe for galaxy formation propounds [see discussion 
in Keres et al. (2005) and references therein]. An important 
point is that our simulations indicate that this formation 
scenario follows from simple physical principles in the con- 
text of the current ACDM scenario. 



5 THREE DIMENSIONAL STRUCTURE OF 
THE BARYONIC CONSTITUENT 

A quantitative description of ELO mass distributions is 
given by their 3D density profile and the structure their 
constituent particles. We first address the structure of the 
baryonic particles. 

5.1 Three dimensional structure for gas particles 

The gas structure is drawn in Figure 1 for the second more 
massive object formed in a SF-B type simulation. The 3D 
density at a given distance, r, from the centre of the object 
has been calculated by binning on concentric spherical shells 
around r. In this Figure, the line is the density profile of 
dark matter around the object, multiplied by Ob/fim- Points 
represent gas density at the positions of SPH particles, and 
colours stand for gas particle temperatures according with 
the scale at the bottom of the Figure. 

We see in this Figure that very few gas is left at posi- 
tions with r < 30 kpc where stars dominate the mass density, 
that cold gas at r > 30 kpc is dense and clumpy, while hot 
gas (that is, gaseous particles with T > 3 x 10 4 K) is diffuse 



with an almost isothermal component at 100 kpc < r < 400 
kpc, and a warm component at the outskirts of the config- 
uration, reaching outside the virial radius (395.0 kpc). Two 
scales stand out in this configuration: the ELO scale or stel- 
lar component, with a size in this case of ~ 30 kpc, and the 
halo scale, a halo of dark matter of 395.0 kpc. Cold dense gas 
particles are associated in most cases with small dark mat- 
ter haloes (not seen in the Figure) ; both gaseous particles in 
cold clumps and dark matter particles in their (sub) haloes 
are shocked particles, using the terminology of the adhesion 
model (see, for example, Vergassola et al. 1994). The con- 
figuration illustrated by this Figure is generic for ELOs: we 
can distinguish an ELO scale, with typical sizes of no more 
than ~ 10 - 40 kpc, and the halo scale, a halo of dark matter 
typically ten times larger in size. 

5.2 Stellar and gaseous particle orbits 

ELO constituent particles of different kinds travel on orbits 
that have different characteristics. To analyse this point, in 
the upper panel of Figure 2 we plot, for each star particle 
and each gaseous particle of a typical ELO, the cosine of the 
angle formed by its position (fl) and its velocity (vl) as a 
function of n . Positions and velocities have been taken with 
respect to the centre of mass of the main baryonic object. 
In this plot radial orbits have cosines = ±1, while circular 
orbits have cosines = 0. Starred (circular) symbols stand for 
stellar (gaseous) particles. We see that cold gas particles at 
r < 4 kpc form a disc in coherent circular motion; filled 
(open) symbols represent particles in corotation (counter- 
rotation) with respect to this small disc. We can also see 
that stellar particle orbits at < 3 kpc scales do not show 
any preference, while those further away, as well as gaseous 
particles outside the disc, show a slight tendency to be on 
radial orbits providing anisotropy to the velocity dispersion. 
Stellar particles constitute a disordered or dynamically hot 
component, showing an important velocity dispersion, and, 
also, in some cases, a coherent net rotation. In §7 these issues 
will be addressed in detail. 



6 THE DARK MATTER AND BARYONIC 
MASS DISTRIBUTIONS 

6.1 Dark matter profiles 

Spherically averaged dark matter density profiles of relaxed 
haloes formed in N-body simulations have been found to be 
well fitted by analytical expressions such that, once rescaled, 
give essentially a unique mass density profile i.e., a two pa- 
rameter family. These two parameters are usually taken to 
be the total mass, M v i r , and the concentration, c or the en- 
ergy content, E. These two parameters are, on their turn, 
correlated (i.e., the mass-concentration relation, see, for ex- 
ample, Bullock et al. 2001; Wechsler et al. 2002; Manrique et 
al. 2003) because the assembly process implies a given cor- 
relation between M v j r and E. Different authors propound 
slightly different fitting formulae, see Einasto (1965, 1968, 
1969) or Navarro et al. 2004, Hernquist 1990 (Hern90), 
Navarro, Frenk & White 1995, 1996 (NFW), Tissera & 
Dommguez-Tenreiro 1998 (TD), Moore et al. 1999 and Jing 
& Suto 2000 (JS), that can be written as: 
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Figure 1. 3D gas (points) and dark matter (blue line) density for a typical ELO. Note the dense cold gas clumps embedded in the 
diffuse hot gas component. See text for an explanation. 



dark/ \ aver c Pfr/^h) 

Ph (r) = Pa x Q ,\ 
35(c) 



(2) 



where PA Cr is the average density within the virial radii, 
c = r v ir/a,h is the so-called concentration parameter, and 



p(y) = y a (i + y) 8 , 



(3) 



where (a,/3) = (1, 3) for Hern90; (a, (3) = (1, 2) for 
NFW; (a,/?) = (2, 2) for TD, and (3 = 3 - a, with a left 
free, for the general formula found by Jing & Suto (2000) 
(note that NFW can be considered as JS with a = 1). In 
these fitting formulae a is the inner slope (r « ah), the 
outer slope (r » ah) is a + (3 (3 for JS or NFW), so that 
ah characterises the scale where the slope changes. Other 
interesting scale is r_2, the r value where the logarithmic 
slope, din p/dlnr = —2. We have r_2 = ah(2 — a)/(a+/3— 2) 
for a profile given by Eq. 3, with r_2 = a/i(2 — a) for JS 
and r_2 = for NFW. Navarro et al. (2004) propound a 
different fitting formula of the form: 



p(y) =e X p(-2 AtJ / 1 ^). 



(4) 



where dlnp/dlnr = —2(r/ah) 1 ^ and r_2 = a^. Note 
that this last fitting formula is similar to the Sersic formula 
(Eq. 1), as Merritt et al. (2005) first pointed out. It was first 
used by Einasto (1965, 1968, 1969), see also Einasto & Haud 
(1989), so that we will refer to it as Einasto model (Eina), 



in consistency with the terminology used by other authors 
(Merritt et al. 2006). 

The g(c) functions can be written as: 



g(y) = y 2 /2(y + l) 2 (HernQO) 

g(y) = ln(y + l)-y/(y + l) (NFW) 

g(y) = 9y/(l + y) (TD) 



(5) 
(6) 
(7) 



g(y) = (3 - ary- Q 2 Fi(3 - a, 3 - a, 4 - a, -y) (JS) (8) 

(9) 



g(y) = i(2 M ) 1 ~ 3 ' J 7(3 M ,2 / V/'<) (Ezna) 



where 2^1 is the hypergeometric function and 7 is the 
lower incomplete gamma function. 

When processes other than gravitational are involved in 
mass assembly (for example, cooling or heating), the dark 
matter density profiles could be modified (see Blumenthal 
et al. 1986; Dalcanton et al. 1997; Tissera & Dommguez- 
Tenreiro 1998; Gnedin et al. 2004). To analyse this point, in 
Figure 3 we plot the dark matter density profiles for several 
typical ELOs, along with their best fit to different analytical 
profiles. The optimal fit has been obtained by minimising the 
statistics: 



X 2 = E£ 1 [logM(< n) - logJWSg(< n)] 2 /N 



(10) 
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Figure 2. Upper panel: the cosine of the angle formed by the position and the velocity vectors for each gaseous (green circles) and 
stellar (starred red symbols) particle belonging to a typical ELO. Filled (open) symbols stand for particles in (counter) corotation with 
the small inner disc 



where Mblo(< r * s the ELO dark matter mass within 
a sphere of radius r% centred at its centre of mass, M(< 
n) is the integrated mass density profile corresponding to 
the different formulae above, and the virial radii r v i r have 
been taken as outer boundaries of the fitting range. The 
integrated dark matter density profiles have been used as 
fitting formulae instead of the dark matter density profiles 
themselves because these latter are binning dependent. An 
updated version of the MINUIT software from the CERN 
library has been used to make these fits as well as the others 
in this paper. 

Note in Figure 3 that the quality of the fits differs from 
one analytical profile to another. To quantify this effect, in 
Figure 4 we plot the distributions of the \ 2 P er d.o.f. statis- 
tics, normalised to (logM v j r ) 2 , resulting from the fits to the 
different profiles above, except for Hern90 one whose results 
are generally poorer. We see that the lower x 2 per d.o.f. val- 
ues generally correspond to either the Eina or the JS profiles, 
with the TD profiles in the third position. In Figure 5 we 
draw the values of the fi (for Eina profiles) and a (for JS 
profiles) slopes corresponding to the optimal fits of SF-A 
sample DM haloes. A slight mass effect can be appreciated 
with lower mass ELOs having steeper DM haloes than more 



massive ones, presumably due to a more important pulling 
in of baryons onto dark matter as they fall to the ELO centre 
with decreasing ELO mass. That is, massive haloes are less 
concentrated than lighter ones, i.e., the mass-concentration 
relation. In any case, the profiles are always steeper than 
a = 1 (i.e., the NFW profile; see Mamon & Lockas 2005a). 

To further analyse this effect, we plot in Figure 6 the 
P-2 density parameter versus the r_2 scale obtained from 
fits to the Einasto model. Blue triangles are measurements 
by Navarro et al. (2004) onto haloes formed in N-body sim- 
ulations and the green line is their best fit. We see that at 
given r_2, p-2 is higher in our hydrodynamical simulations 
than in those of Navarro et al. (2004), presumably due to the 
pulling in of dark matter by baryon infall. We also see that 
at given M v i r , r_2 is shorter in hydrodynamical simulations 
than in purely gravitatory ones, by the same reason. 



6.2 Projected stellar mass density profiles 

Authors now agree that the Sersic law given in Eq. 1 (Sersic 
1968) is an adequate empirical representation of the optical 
surface brightness profiles of most ellipticals (see, for exam- 
ple, Caon et al. 1993; Bertin et al. 2002). Assuming that the 
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Figure 3. Dark matter density profiles (black full line) for several typical ELOs from SF-A and SF-B samples along with their best 
fits to different analytical profiles: NFW (red point line), TD (blue long-dashed line), JS (green short-dashed line) and Eina (magenta 
point-dashed line). 



stellar mass-to-light ratio i star does not appreciably change 

with ELO projected radius R §, the projected stellar mass 
profile, E star (i?) can be taken as a measure of the surface 
brightness profile and be written as 

E star (#) = 7 star / light ( J R). (11) 

One can then expect that E star (7?) can be fitted by a 
Sersic-like law. This is in fact the case as shown in Fig- 
ure 7 for several typical ELOs drawn from both SF-A and 
SF-B samples (see Kawata & Gibson, 2005, for a similar re- 
sult concerning one virtual elliptical galaxy). Some remarks 



Hereafter we will use capital R to mean projected radii 



on how our fits have been made are in order. First, the 
E star (7?) p ro gi es nave been calculated by averaging on con- 
centric rings centred at the projection of the centre-of-mass 
of the corresponding ELO. Three projections along orthogo- 
nal directions have been considered for each ELO. Also, be- 
cause these projected densities are binning dependent and 
somewhat noisy, the integrated projected mass density in 
concentric cylinders of radius R and mass 



M%f(R) = 2tt / E star (R')R' dR' 



(12) 



has been used as a fitting function, instead of E star (J?) 
itself. Concerning the fitting range, we have adopted an 
outer boundary i? m ax such that the corresponding surface 
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Figure 4. The distributions of the x 2 per d.o.f., normalised to 
the logarithm of their respective mass square, for the fits of the 
DM density profiles of ELO haloes (SF-A and SF-B samples) to 
different analytical profiles 



brightness l light (7? ma x) (see Eq. 11) gives the standard value 
of \ib = 27 mag arcsec -2 . The values for the stellar mass-to- 
blue-light 7^ ar span a range from -f^ = 2 to 12, depending 
on the details of its determination (see discussion in Mamon 
& Lockas 2005a), and best values of 7s tar = 5 to 8. Their 
geometric mean 7fj ar = 6.3 has been used to make the fits 
drawn in Figure 7, but the results of the fit do not signifi- 
cantly depend on the particular 7# ar value used within its 
range of best values. 

An interesting result is that the values of the shape 
parameter n we have obtained are consistent with obser- 
vations, including their correlations with the effective radii 
i?e ght , luminosity L and velocity dispersion (Caon et al. 
1993; Prugniel & Simien 1997; Graham 1998; Marquez et 
al. 2000; D'Onofrio 2001; Trujillo et al. 2001; Vazdekis, Tru- 
jillo & Yamada 2004; Graham et al. 2006), as illustrated in 
Figure 8. In this Figure we plot the shape parameter n versus 
the ELO projected stellar half-mass radii, R% defined by 
the condition that M c s * ar (Rf£ 



encloses half the total stel- 
lar mass of the system; assuming that 7§ ar does not depend 
on R, we will have -Rjfbo — -Re ght - Blue triangles are data on 
n and i?L ight from D'Onofrio (2001). Note that a slight effect 
resulting from the different SF parametrization in SF-A and 
SF-B sample ELOs is apparent in this plot, mainly due to 
the smaller sizes of SF-B sample ELOs as compared with 
their SF-A sample counterparts. 



6.3 Baryonic three-dimensional mass density 
profiles 

In the last section it has been shown that the projected 
stellar mass density profiles are adequately described by the 
standard Sersic profiles, that is, that they are consistent with 
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Figure 5. Left panel: the optimal inner slope a of the general Jing 
& Suto profile for the DM haloes of ELOs (green filled squares) 
and the /i coefficient of the Einasto analytical profile (magenta 
filled circles), versus their virial mass for SF-A sample ELOs. 
Right panel: zoom of the a versus virial mass plot to clarify the 
mass effect. 



observational data. We now analyse the three-dimensional 
mass density profiles of baryons. 

We first analyse the baryon distribution at the ELO 
scale, where the main contribution to the mass density 
comes from stars. We lack of any observational input on how 
the three-dimensional stellar-mass density profiles p star (r) 
can be, except for a deprojection of the Sersic profiles (Prug- 
niel & Simien 1997; Lima Neto, Gerbal & Marquez 1999). 
In Figure 9 we plot /5 star (r) for ELOs in the SF-A sample. 
Different colours have been used for ELOs in different mass 
intervals and a clear mass effect can be appreciated in this 
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Figure 6. The p_2 density parameter versus the r_2 scale ob- 
tained from fits to the Einasto model, for ELOs in both the SF-A 
sample (filled black circles) and the SF-B sample (open red cir- 
cles). Blue triangles are measurements by Navarro et al. (2004), 
onto haloes formed in N-body simulations, with its fit by Mamon 
& Lockas 2005a (green line). Numbers correspond to the loga- 
rithms of the virial masses (in units of Mq) of haloes formed in 
different simulations, according with their respective colours 
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Figure 7. Projected stellar mass density profiles for different 
ELOs along with their best fit by a Sersic law. The correspond- 
ing shape parameter best values and minimal x 2 per-degree-of 
freedom are also shown 
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Figure 8. The Sersic shape parameter n versus the projected 
stellar half-mass radii for SF-A sample (black filled circles) and 
SF-B sample (red open circles) ELOs. For each ELO, the results 
of projections along three orthogonal directions are shown. Blue 
filled triangles are data on n and Re from D'Onofrio (2001) 



[H] 





Figure 9. Three-dimensional stellar mass profiles for ELOs in the 
SF-A sample: full green lines, ELOs with M vir < 1.5 X 10 12 M ; 
orange point-dashed lines, ELOs with 1.5 X 10 12 M Q < M vir < 
5 X 10 12 M Q ; blue dashed lines: ELOs with M vil > 5 X 10 12 M . 
The stellar mass density profiles show homology breaking 
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Figure 11. The \ 2 statistics (Eq. 10) corresponding to the fits of 
the stellar mass profiles for ELOs in the SF-A sample Einasto pro- 
files (magenta filled circles) and JS profiles (green filled squares) 



Figure, and particularly so at the inner regions, where at 
fixed r/r v i r the stellar-mass density of less massive ELOs 
can be a factor of two or so higher than that of more mas- 
sive ones. This means that the mass homology is broken in 
the three-dimensional stellar mass distribution. 

To quantify the stellar three-dimensional mass density 
profiles of ELOs, they have been fit to JS and Einasto an- 
alytical formulae through the statistics defined in Eq. 10 
where Mg"o(< r») has been replaced by the ELO stellar 
mass within a sphere of radius r». The quality of the fits is 
illustrated in Figure 10, and in Figure 11 the values of the 
X 2 p.d.o.f. statistics are given, normalised to logM^f 1 . Both 
Figures show that these profiles describe adequately well 
the spherically-averaged stellar mass distribution in three 
dimensions, even if with very small r_2 values. 

To study the possibility that the homology in the dark- 
versus bright-mass distribution is also broken, the stellar-to- 
dark density ratio profiles 

/p tar M = p star M/p dark M (13) 

are plot versus either the radii (Figure 12 upper panel) 
or the radii normalised to virial radii (Figure 12 lower panel). 
We see that there is, in any case, a clear mass effects at the 
inner regions, with the stellar mass distribution relative to 
the dark mass one less concentrated with increasing ELO 
mass. For example, in Figure 12 we see that the fraction of 
the ELO virial volume where /p tar (r) > 1, is smaller as the 
ELO mass grows; also, at fixed r/r v i r , / star (r) increases with 
decreasing ELO mass. So, the homology is broken in the 
three-dimensional stellar-to-dark mass distribution, a fact 
that could be important to explain the tilt of the observed 
FP (see Onorbe et al. 2005, 2006). 

To further analyse this point and make the comparison 
with observational results easier, the dark-to-stellar mass ra- 
tio profiles, M dark (< r)/M tot (< r), are drawn in Figure 13 
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Figure 12. The stellar-to-dark mass profiles for ELOs in the SF- 
A sample. In the lower panel, radii are normalised to the virial 
radii. Line types and colours are the same as in Figure 9 



for the same ELOs, with the radii in units of the three di- 
mensional stellar half-mass radii T We see that there is, in 
any case, a positive gradient, and again a clear mass effect, 
with a tendency of the dark matter fraction at fixed val- 
ues of r/r^bo t° be higher as the mass scale increases. To 
be more quantitative and compare with observational data, 
we plot in Figure 14, upper panel, the fraction of dark-to- 
total masses at r/R B ^ T = 1 for ELOs in the SF-A and SF-B 
samples, versus their stellar masses. The differences among 
results for both samples come from the smaller i?c!bo values 

^ The effective or stellar half-mass radii r**^ are defined as the 
radii of the spheres enclosing half the ELO stellar mass. This is 
the relevant three-dimensional size parameter at ELO scales. See 
Onorbe et al. (2006) for details 
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Figure 10. The stellar mass profiles for 4 typical ELOs in SF-A sample (black continuous lines) and their optimal fits to Einasto profiles 
(magenta point-dashed lines) and JS profiles (green dashed lines) 



of SF-B sample ELOs as compared with their SF-A counter- 
parts (see details in Onorbe et al. 2006). Blue triangles with 
error bars are results from integral field SAURON data and 
models by Cappellari et al. 2006. We see that these empiri- 
cal determinations of the dark matter fraction at the centre 
of ellipticals is consistent with the values found in ELOs of 
both samples, including its growth with the mass scale. 

In the lower panel of Figure 14 we give the gradients 
of the Af dark (< r)/M star (< r) profiles as a function of their 
stellar masses. Blue triangles with error bars are the empir- 
ical mass-to-light gradients as determined by Napolitano et 
al. 2005 for EGs with isophotal shape x 100 < 0.1, that 
is, boxy ellipticals. We have used as inner and outer radii 

/ star, in a r l / star.out * l 1 j_i l 

r/r c bo = 0.5 and r/r c bo = 4, roughly the average val- 
ues of the inner and outer radii these authors give in their 



Table 1. We see that there is a mass effect and that our 
results are consistent with those found by these authors in 
the range of stellar mass values our samples span, especially 
when we consider that ELOs in our samples are boxy (see 
§7.2). A SF effect in the stellar mass distribution also ap- 
pears in Figure 14, again due to the compactness of the SF-B 
sample ELOs relative to their SF-A sample counterparts. 

We now turn to analyse the baryon space distribution 
at halo scales. To have an insight on how baryons of any 
kind are distributed relative to the dark matter at the halo 
scale and beyond, the baryon fraction profile 

/ ab (r)=P ab M/p t0t M, (W) 

where "ab" stands for baryons of any kind (i.e., stars, 
cold gas and hot gas) and "tot" stands for matter of any 
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Figure 13. The fraction of dark-to-total mass profiles, Af dark (< 
r)/M tot (< r) for ELOs in the SF-A sample; different colours 
stand for different mass ranges, as in Figure 9; radii are normalised 
to the 3D stellar half-mass radii 



kind (i.e., dark plus baryons of any kind), is drawn in Fig- 
ure 15 for ELOs in the SF-A sample (black full lines) and 
in the SF-B sample (red point lines) in the same range of 
virial mass, 1.5 x 10 12 M Q < M vil < 5 x 1O 12 M . Despite 
individual characteristics, the / ab (r) curves show a typical 
pattern in which their values are high at the centre, then 
they decrease and have a minimum lower than the global 
value, /cosmo = Qb/Qm = 0.171, at a radius r^ in , then they 
increase again, reach a maximum value and then they de- 
crease and fall to the /^ smo value at a rather large r value, 
larger than the corresponding virial radii. This result, i.e., 
that EGs are not baryonically closed, is also indicated by 
recent X-ray observations (Humphrey et al. 2006). Notice 
(Figure 16) that the increase of / ab (r) at r > r ab n is mainly 
contributed by hot gas, almost absent at r < r^ in , indicat- 
ing that rjfij, separates the (inner) region where gas cooling 
has been possible from the (outer) region where gas has not 
had time enough to cool in the ELO lifetime. Note also in 
Figure 16 that an important amount of hot gas is outside 
the spheres of radii r v i r , that is, it is not bound to the self- 
gravitating configuration defined by the ELO halo. In fact, 
the mass of hot gas increases monotonically up to r ~ 4r v ir, 
and maybe also beyond this value, but it is difficult at these 
large radii to properly dilucidate whether or not a given hot 
gas mass element belongs to a given ELO or to another close 
one (to alleviate this difficulty, only those ELOs not having 
massive neighbours within radii of 6xr v i r have been consid- 
ered to draw this Figure) . Another important result is that 
the hot gas mass fraction, relative to the cold mass fraction 
at the ELO scale, increases with M v i r at given r/r v i r , and 
the differences between massive and less massive ELOs can 
be as high as a factor of ~ 2 at r/r v i r < 4. We see that, in 
massive ELOs, this excess of baryons in the form of hot gas 
at the outer parts of their configurations, compensates for 
the lack of baryons in the form of stars at the ELO scales. 
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Figure 14. Upper panel: the fraction of dark-to-total mass at 
r/il st ? r = 1 versus the ELO stellar masses. Filled (open) sym- 
bols: ELOs in SF-A (SF-B) sample. Points with error bars are the 
values corresponding to the SAURON sample of ellipticals. Lower 
panel: the gradients of the M dark (< r)/M star (< r) profil es as a 
function of their stellar masses; blue triangles with error bars are 
the empirical mass-to-light gradients as determined by Napoli- 
tano et al. 2005 for galaxies with the (14 X 100 shape parameter 
lower than 0.1 (that is, for boxy ellipticals) 



6.4 Total three-dimensional mass density profiles 

We now address the issue of the total mass (i.e., baryonic 
plus dark) density profiles. In Figure 17 they are drawn 
for ELOs in the SF-A sample (upper panel) as well as for 
those in the SF-B sample (lower panel). In both cases, the 
profiles corresponding to ELOs in different mass intervals 
have been drawn with different line and colour codes. Some 
important results are that i), they are well fit by power- 
law expressions p tot (r) oc r~ 7 in a range of r/r v i r values 
larger than two decades, ii), the slope of the power-law in- 
creases with decreasing ELO mass, and, iii) a slight SF ef- 
fects appears, but only at the very inner regions, with SF- 
B sample ELOs showing a worse fit to a power law than 
their SF-A counterparts. Koopmans et al. 2006 have also 
found that the total mass density profiles of their massive 
(< o" ap >= 263 km s _1 ) lens EGs can be fit by power-law ex- 
pressions within their Einstein radii (< i^Einst >= 4.2 ± 0.4 
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Figure 15. Baryon fraction profiles for ELOs in SF-A sample 
(black full lines) and SF-B sample (red point lines), in the same 
range of virial mass, 1.5 X 10 12 M Q < A/ vir < 5 X 10 12 M Q 
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Figure 16. The Af h s(< r)/M^ profiles for typical ELOs. 
M h s(< r) is the mass of hot gas within a sphere of radius r. 
Orange point-dashed lines: ELOs with 1.5 X 10 12 < Af v i r < 
5 X 1O 12 M ; green full lines: ELOs with M vir < 1.5 X 1O 12 M . 
Only isolated ELOs have been considered. 



kpc, with < R E inst/Re Sht >= 0.52 ± 0.04, i.e., the inner re- 
gion), whose average slope is < 7 >= 2.0llo.o3 ± 0.05 (68 
percent C.L.), with an intrinsic scatter of 0.12. To elucidate 
how well the total mass density profiles of ELOs compare 
with these results, in Figure 18 we plot the slopes 7 for 
ELOs, as well as for SLACS lens ellipticals (Table 1, Koop- 
mans et al. 2006), versus their central L.O.S. stellar velocity 
dispersions. The fitting range for ELOs used to draw this 
Figure is r < Rl%. We see that results for ELO and SLACS 
lens galaxy samples are consistent in the range of velocity 
dispersion values where they coincide. 



7 KINEMATICS 

7.1 Three-dimensional velocity distributions 

Shapes and mass density profiles (i.e., positions) are related 
to the 3D velocity distributions of relaxed EGs through the 
Jeans equation (see Binney & Tremaine 1987). Observation- 
ally, informations on such 3D distributions are not available 
for external galaxies, only the line-of-sight velocity distribu- 
tions (LOSVD) can be inferred from their spectra (Binney 
& Tremaine 1987; van der Marel & Franx 1993; de Zeeuw 
& Franx 1991). The complete six dimensional phase space 
informations for each of the particles sampling the ELOs 
provided by numerical simulations, allow us to calculate the 
3D profiles for the velocity dispersion, o"3d(V), as well as the 
circular velocity profiles, V C i T (r) . In Figure 19 we draw the 
Vri r (r) profiles (full line), as well as their dark matter (short- 
dashed line) and baryonic contributions (stars, long-dashed 
line; stars plus cold gas, point line). 

The Vcir(r) profiles provide another measure of ELO 
mass distribution. We note in Figure 19 that the baryon 
mass distribution is more concentrated than the dark matter 
one due to energy losses by the gaseous component before 



being transformed into stars. This is a general property of 
the circular velocity profiles of the ELO samples. Moreover, 
objects in SF-B sample are more concentrated than their SF- 
A sample counterparts, because of the SF implementation: 
the amount of baryons at their central volumes relative to 
dark matter is always lower in SF-A than in SF-B objects; 
this is a small scale effect as r ~ 30 kpc or r ~ 40 kpc radii 
enclose roughly similar amounts of baryons or dark matter 
in any cases. 

In Figure 20 we draw, for the same ELOs, the <T3d(»") 
profiles as measured by stars, (f|^ r (r), starred symbols and 
short-dashed lines), and by dark matter, (crip (r), open cir- 
cles and long-dashed lines) as proof particles in the overall 
potential well. These profiles are in any case decreasing out- 
wards, both for the dark matter and for the stellar com- 
ponents. An outstanding result illustrated by Figure 20 is 
that &3 D (r) is always higher than cr^ T (r) (because stars 
are made out of cooled gas), with (jfn{r)/af^(r) ~ 0.8, in 
consistency with the values found by Loewenstein (2000) on 
theoretical grounds and by Dekel et al. (2005) from pre- 
prepared simulations of mergers of disc galaxies. This is 
the so-called kinematical segregation (Saiz 2003, Saiz et al. 
2004). To further analyse this issue, in Figure 21 we plot 
the af3 r (r)/a?S k (r) ratios for the ELOS in both the SF- 
A sample and in the SF-B sample, with different colour 
and line codes depending on the ELO mass range. We see 
that the kinematical segregation does not show either a clear 
mass dependence, or a radial dependence. Moreover, the SF 
parametrization effect is only mild. 

Another relevant quantity is the anisotropy of the 3D 
velocity distributions of the ELO sample, defined as: 

/3ani = l-^, (15) 

where a T and a t are the radial and tangential veloc- 
ity dispersions (<rj? = <r 2 + a 2 ,), relative to the centre of the 
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Figure 17. The total mass density profiles for ELOs in the SF-A 
sample (up) and in the SF-B sample (down). Green full lines: 
ELOs with M v j r < 1.5 X 1O 12 M0; orange point-dashed lines: 
ELOs with 1.5 X 1O 12 M < Af vir < 5 X 10 12 M ; blue dashed 
lines: ELOs with 5 X 1O 12 M < A/ vir . The violet long-dashed 
lines are the one sigma interval for the slope resulting from fits to 
power-law profiles of lens ellipticals from Koopmans et al. (2006). 



object. The anisotropy profile, /3 an i(?"), is represented in Fig- 
ure 20 for typical ELOs in the sample, for their dark matter 
and stellar particle components. The anisotropy is always 
positive (i.e., an excess of dispersion in radial motions), the 
profiles are almost constant, except at the innermost regions, 
and the stellar component is always more anisotropic than 
the dark matter one, presumably as a consequence of the 
mergers involved in the ELO mass assembly (see §4 and 
DTal06). In fact, the characteristics of the stellar anisotropy 
profiles (roughly constant and /3anT( r ) — 0.5 in most cases) 
are consistent with those found by Dekel et al. (2005), where 
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Figure 18. The logarithmic slopes corresponding to the total 
mass density profiles for ELOs in the SF-A (black filled circles) 
and the SF-B samples (red open circles), versus their central 
L.O.S. stellar velocity dispersions. Triangles with error bars cor- 
respond to data on SLACS lens ellipticals, as given in Table 1 of 
Koopmans et al. (2006) 



they conclude that large radial anisotropy is generic to the 
stellar component of merger remnants of any kind. 



7.2 Stellar LOS velocity and velocity dispersion 
profiles 

The two last Figures provide an illustration of the general 
characteristics of the lower-order moments of the 3D veloc- 
ity distribution. The profiles plot in these Figures are not 
observationally available, but only the lower-order moments 
of the LOSVD are. We have measured the stellar line-of- 
sight velocity and the stellar velocity dispersion profiles, 
Vios T (R) an d crf^ T (R), along three orthogonal projections 
for all ELOs in the sample. To mimic observational tech- 
niques used in stellar kinematics of EGs, we have measured 
these profiles along the major and minor axes of projected 
ELOs, where the major axis is defined as that orthogonal to 
the ELO spin vector projected on the plane normal to the 
LOS, and the minor axis is parallel to the spin projection. 
We have found that in some cases ELOs do indeed show a 
clear rotation curve, while in most cases the rotation is only 
modest or even very low, as illustrated in Figure 22 and in 
Figure 23, respectively. 

To quantify the amount of rotation in ELOs and its 
possible dependence on the mass scale, in Figure 24 we plot 
the ratios c rot = Vmin/ (V^j + V^ in ) 1//2 as a function of the 
ELO virial masses, for ELOs in both the SF-A and the SF- 
B samples (Vmaj and Vmin are the maximum values of the 
Viog(i?) profile when measured along the major and the 
minor axes, respectively). When the ELO shows a clear ro- 
tation curve, V m in is much lower than Vmaj , and the c rot ratio 
is low; by contrast, when the rotation is unimportant, then 



© 0000 RAS, MNRAS 000, 000-000 



Bright and Dark Matter in Ellipticals 17 




0716.27875.850.534.303 8722.28648.228.828.311 

Figure 19. The circular velocities profiles of two typical ELOs in the SF-A sample (upper panels) and their SF-B sample counterparts 
(lower panels). Black continuous line: total mass; blue short-dashed line: dark matter contribution; green long-dashed line: stellar mass 
contribution; red point line: cold baryon contribution 



Vmin o; Vma.j and Crot ~ 0.7. For a given ELO, the c ro t value 
depends on the direction taken as LOS direction, in such a 
way that it is maximum when the ELO spin is taken as LOS 
direction, and minimum when the LOS direction is normal 
to the ELO spin vector, that is, when rotation stands out. 
This is the LOS direction taken to draw this Figure, where 
we see that there is not a clear mass dependence, that most 
ELOs are in between the two situations described above 
and that the values of the c rot ratio of ELOs are typical 
of boxy ellipticals (see, for example, Binney & Merrifield, 
figure 4.39). 

We now comment on the major axis LOS stellar ve- 
locity dispersion profiles of ELOs (Figure 22 and in Fig- 
ure 23). Their most outstanding feature is the decrease of 
the ofo^ r (i?) profiles in some cases and particularly so along 
some LOS directions at large R. These profiles are suited 
to compare with stellar kinematics data. In other cases, for 
example to compare with planetary nebulae data, the LOS 
velocity dispersion profiles must be calculated by averaging 
over the LOS velocities of stars placed within cylindrical 

© 0000 RAS, MNRAS 000, 000-000 



shells, with their axes in the LOS direction. Figure 25 is a 
plot of such profiles normalised to crf^ 1 (Rl^bo) for the SF-A 
sample ELOs; each panel corresponds to a different orthog- 
onal projection. 

To make clearer the decline of the <xf*g r (i2) profiles, in 
Figure 26 we plot, at different R values, the averages of the 
stacked profiles shown in Figure 25 with their dispersions 
(green points and error bars), as well as the averages of the 
profiles corresponding to young stars (age < 3 Gyears, or- 
ange squares and error bars), normalised for each ELO to 
their corresponding crf*"^*^.). The decline of these ve- 
locity dispersion profiles can be clearly appreciated, as well 
as the slightly larger decline of the profiles corresponding 
to the younger stellar populations. These results are con- 
sistent, within their dispersions, with that shown by Dekel 
et al. (2005) in their figure 2 (lower panel). They are also 
marginally consistent with the decline shown by PN data in 
the NGC 821, NGC 3379, NGC 4494 and NGC 4697 galaxies 
(Romanowsky et al. 2003; Romanowsky 2006). Note, how- 
ever, that our ELOs are boxy, while the 04 x 100/a shape 
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Figure 20. The (T3Y>(r) profiles of two typical ELOs in the SF-A sample (upper panels) and their SF-B sample counterparts (lower 
panels). Also shown are the anisotropy profiles /3 a ni(0- Long-dashed lines: dark matter; short-dashed lines: stars 
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Figure 21. The o"| 1 g r (r)/(Tg| ) rk (r) ratio profiles for ELOs in SF-A (upper panel) and SF-B (lower panel) samples. Different colour and 
line codes stand for ELO mass intervals, as in Figure 9 
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Figure 22. Upper panel, full line: the major axis stellar LOS 
velocity profile along the spin direction for an ELO in SF-A sam- 
ple. Point and dashed lines: same as the continuous line taking 
the LOS direction normal to the ELO spin vector. This particular 
ELO rotates. Lower panel: same as the upper panel for the major 
axis LOS stellar velocity dispersion profiles 
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Figure 23. Same as the previous Figure for another ELO. In this 
case, the rotation is only mild. 



parameters for these galaxies are 2.5, 0.2, 0.3 and 1.4, re- 
spectively, that is, they are rather discy ellipticals. 



8 SUMMARY, DISCUSSION AND 
CONCLUSIONS 

8.1 Summary: method and main results 

This paper belongs to a series aimed at studying galaxy 
formation in a cosmological context through hydrodynam- 
ical simulations. Here we present an analysis of the sam- 
ple of elliptical-like-objects (ELOs) formed in ten different 
cosmological simulations, run within the same global flat A 
cosmological model, roughly consistent with observations. 
The normalisation parameter has been taken slightly high, 
as = 1.18, as compared with the average fluctuations of 
2dFGRS or SDSS galaxies, to mimic an active region of the 
Universe. Newton laws and hydrodynamical equations have 
been integrated in this context, with a standard cooling al- 
gorithm and a SF parameterisation through a Kennicutt- 
Schmidt-like law, containing our ignorance about its details 
at sub-kpc scales, and where subresolution processes affect- 
ing SF are implicitly taken into account through the values 
given to these parameters. No further hypotheses to model 
the assembly processes have been made. Individual galaxy- 
like objects naturally appear as an output of the simulations, 
so that the physical processes underlying mass assembly can 
be studied. Five out of the ten simulations (the SF-A type 
simulations) share the SF parameters and differ in the seed 
used to build up the initial conditions. To test the role of 
SF parameterisation, the same initial conditions have been 
run with different SF parameters making SF more difficult, 
contributing another set of five simulations (the SF-B type 
simulations). ELOs have been identified in the simulations 
as those galaxy-like objects at z — having a prominent, dy- 
namically relaxed spheroidal component made out of stars, 
with no extended discs and very low gas content. These stel- 
lar component is embedded in a dark matter halo that con- 
tributes an important fraction of the mass at distances from 
the ELO centre larger than ~ 10 — 15 kpc on average, within 
which some clumps made out of cold dense gas and stars, 
associated in some cases with dark matter substructures, or- 
bit. No ELOs with stellar masses below 3.8 x 10 10 M or 
virial masses below 3.7 x 10 11 M have been found that 
met the selection criteria (see Kauffmann et al. 2003 for a 
similar result in SDSS galaxies and Dekel & Birnboim 2006, 
and Cattaneo et al. 2006 for a possible physical explana- 
tion). ELOs have also an extended halo of hot, diffuse gas. 
Stellar and dark matter particles constitute a dynamically 
hot component with an important velocity dispersion, and, 
except in the very central regions, a positive anisotropy. 

The first step in the program of studying the origins of 
EGs through self-consistent simulations, is to make sure that 
they produce ELO samples that have counterparts in the 
real local Universe. This objective has been partially fulfilled 
in previous works. An analysis of the structural and dynam- 
ical ELO parameters that can be constrained from obser- 
vations has shown that they are consistent with those mea- 
sured in the SDSS DR1 elliptical sample (Saiz et al. 2004), 
including the FP relation (Onorbe et al. 2005; Onorbe et al. 
2006). Also, ELO stellar populations have age distributions 
with the same trends as those inferred from observations, 
i.e., most stars have formed at high z on short timescales, 
and, moreover more massive objects have older means and 
narrower spreads in their stellar age distributions than less 
massive ones (DSS04). 
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Figure 25. LOS velocity dispersion velocity profiles along three different orthogonal projections for ELOs in SF-A sample up to 6 
effective radii. The profiles are normalised to their value at R^bo ^ or eacn ELO. Different colour and line codes stand for ELO mass 
intervals, as in Figure 9 



In this paper we address the important issue of the 
amount and distribution of dark matter in virtual ellipti- 
cals and, in particular, its amount and distribution relative 
to the bright matter distribution. We also address the kine- 
matics of the dark matter component and its relationship 
with the kinematics of the bright matter component. We 
want to answer to the question of the mass and extension 
of dark matter haloes in EGs. To answer to this question, 
we have first compared our virtual results with new obser- 
vational data, obtaining a very satisfactory agreement. To 
be specific: 

• The projected stellar mass profile, E star (7?), can be ad- 
equately fitted by a Sersic-like law. The shape parameter 
values n we have obtained are consistent with observations, 



including their correlations with the ELO luminosity, size 
and velocity dispersion. 

• The fraction of dark-to-total mass inside the projected 
half-mass radii are consistent with the observational ones 
obtained by Cappellari et al. (2006) from SAURON data. 

• The gradients 
of the dark-to-stellar M dark (< r)/M star (< r) profiles as a 
function of their stellar masses, are consistent with those 
observationally found by Napolitano et al. (2005) for boxy 
ellipticals. 

• The total mass (i.e., baryonic plus dark) density profiles 
can be well fit by a power-law expression in a large range 
of r/r v ir values, with power-law slopes that are consistent 
with, within the dispersion, or slightly higher than those 



© 0000 RAS, MNRAS 000, 000-000 



Bright and Dark Matter in Ellipticals 21 



> 
+ 



S 0.1 



10" 



(M ) 



Figure 24. The c ro t ratios as a function of the virial mass for 
ELOs in the samples 



observationally found by Koopmans et al. 2006 for massive 
lens ellipticals within their Einstein radii. 

• The line-of-sight velocity profiles along the major axis 
show, in some cases, a clear rotation, even if in most cases 
the rotation is modest or low. The values of the rotation ratio 
along the major and minor axis (a measure of the rotation 
in ELOs) does not depend on the mass scale 

• The values spanned by the rotation ratios of ELOs are 
typical of boxy ellipticals. 

• The line-of-sight velocity dispersion profiles, ai os (R), 
decline outwards at large R, and the slope slightly increases 
when only the younger stellar populations are considered. 
These profiles are only marginally consistent with data on 
PNs at large radii; but these correspond to discy ellipticals 
while our virtual ellipticals are rather boxy. 

These agreements strongly suggest that the intrinsic 
three-dimensional dark and bright matter mass and velocity 
distributions we get in our simulations might also adequately 
describe real ellipticals. We now summarise our most im- 
portant findings on the three-dimensional mass and velocity 
structure of ELOs: 

• ELOs are embedded in extended massive dark matter 
haloes. 

• The best fits of their spherically-averaged dark matter 
density profiles to usual analytical formulae (Hern90, NFW, 
TD, JS, Eina) are generally provided by the two last for- 
mulae. The quality of the fits is good, so that ELO haloes 
form a two-parameter family where the two parameters are 
correlated. This is consistent with those produced in purely 
N-body simulations. The JS inner slope parameter, a, is al- 
ways higher than the NFW value (a =1). 

• The slope parameters grow as the ELO mass scale de- 
creases, indicating that the halo concentration grows when 
the mass decreases. 

• Halos have suffered from adiabatic contraction. This 



can be made quantitative by comparing the plot of the den- 
sity at the Einasto scale, p-2, versus the scale r_2 = an, 
with the plot provided by Navarro et al. 2004 (results of 
purely N-body simulations). 

• At the ELO scale, most baryons have turned into stars. 
The three dimensional stellar-mass density profiles can be 
fit by Einasto or JS profiles, but with small r_2 values. 

• The mass distribution homology is broken in the stel- 
lar mass as well as in the dark- versus bright-mass distri- 
butions, with the stellar mass distribution relative to dark 
mass one less concentrated with increasing ELO mass. That 
is, massive ELOs miss baryons at short scales as compared 
with less massive ones, when we normalise to the dark mat- 
ter content. This result is consistent with the observational 
ones by Cappellari et al. (2006) from SAURON data, as well 
as by Napolitano et al. (2005) we quoted above. 

• At the halo scale, the baryon fraction profiles have been 
found to show a typical pattern, where their values are 
high at the centre, then they decrease and have a minimum 
roughly at 0.3 < r^ in /r v i T <0.7, well below the global value, 
Qb/^im = 0.171, then they increase again, reach a maximum 
value and then they decrease and fall to the global flt/Qm 
value well beyond the virial radii r v i r . This suggests that the 
baryons that massive ELOs miss at short scales (stars) are 
found at the outskirts of the configuration as diffuse hot gas. 
This result could reflect the presence of a stable virial shock 
that prevents gas infall more efficiently as mass increases 
(Dekel & Birnboim 2006). 

• Concerning kinematics, ELO velocity dispersion pro- 
files in three dimensions are slightly decreasing for increas- 
ing r, both for dark matter and stellar particles, (Tg D (r) and 



• The dark and bright matter components of ELOs are 
kinematically segregated, as we have found that (o"3 D (r)) 2 ~ 
(1.4 - 2) {af£ T (r)) 2 , confirming previous results (Saiz 2003; 
Loewenstein 2000; Dekel et al. 2005). This is so because stars 
are formed from gas that had lost energy by cooling. 

• This kinematical segregation does not show any clear 
mass or radial dependence. 

• The anisotropy is always positive (i.e., an excess of 
radial motions) and almost non-varying with r inside the 
ELOs. Recall, however, that ELOs have been identified as 
dynamically relaxed objects: there are not recent mergers in 
our samples. 

• The stellar component generally shows more anisotropy 
than the dark component, maybe derived from the radial 
motion of the gas particles that gave rise to the stars. 



8.2 Possible resolution effects 

To make sure that the results we report in this paper are 
not unstable under resolution changes, a control simulation 
with 128 3 dark matter and 128 3 baryonic particles, a gravi- 
tational softening of e — 1.15 kpc and the other parameters 
as in SF-A type simulations (the S128 simulation), has been 
run. The results of its analysis have been compared with 
those of a 2 x64 s simulation (the S64 simulation), whose 
initial conditions have been built up by randomly choosing 
1 out of 8 particles in the S128 initial conditions, so that 
every ELO in S128 has a counterpart in the lower resolution 
simulation and conversely. Due to the very high CPU time 
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Figure 26. The SF-A sample average LOS velocity dispersion 
profiles normalised to their values at fijj'bo f° r eacn ELO (green 
points) along with their 1 a dispersions. Orange points and er- 
ror bars: the same for the young stellar particles, with the same 
normalisation 




r[Kpc] 



Figure 27. Three-dimensional stellar mass density profiles for 2 
ELOs identified at z = 1 in Sf 28 (red lines) and their counterparts 
in S64 (black lines). 



requirements for SI 28, the comparison has been made at 
z — 1. The results of this comparison are very satisfactory, 
as Figure 27 illustrates. 

Otherwise, Figures 20 and 21 show that two-body re- 
laxation effects (typically the most stringent requirement 
for convergence) have not been important at least for r 
larger than ~ 1 kpc. In fact, two-body relaxation ef- 



fects cause energy equipartition. But the the values of 
the (jio(r)/aio (r) ratios we have obtained (~ 0.8) ex- 
clude energy equipartition among dark matter and stellar 
particles in ELOs, because such equipartition would de- 
mand afS z (r)/o-3S rk (r) = [m dark /m star ]°' 5 = 2.194, where 

m dark = 1 2 g x 1O 8 M0 and m star = 3 g 7 x 1O 7 M are the 

mass of dark and stellar particles, respectively. This result is 
important because it shows that two-body relaxation effects 
have played no important role in the gravitational interac- 
tion. 



8.3 The physical processes underlying mass 

homology breaking and their observational 
implications 

One of the most important findings of this paper is the ho- 
mology breaking ELO samples show both in the relative 
content and in the relative distribution of the baryonic and 
the dark mass components. As explained in Onorbe et al. 
(2005, 2006), this has as a consequence the observed tilt of 
the Fundamental Plane relation relative to the virial one. 
Which are the physical processes underlying this breaking 
of homology? According with our simulations, they lie in 
he systematic decrease, with increasing ELO mass, of the 
relative amount of dissipation experienced by the baryonic 
mass component along ELO stellar mass assembly (DTal06, 
Onorbe et al. 2006). This possibility, already suggested by 
Bender, Burstein & Faber (1992), Guzman, Lucey & Bower 
(1993) and Ciotti, Lanzoni & Renzini (1996), was first ad- 
dressed through numerical methods by Bekki (1998). He 
studied elliptical formation through pre-prepared simula- 
tions of dissipative mergers of disc galaxies, where the ra- 
pidity of the SF in mergers is controlled by a free efficiency 
parameter Csf. He shows that the SF rate history of galaxies 
determines the differences in dissipative dynamics, so that 
to explain the lack of homology in EGs he needs to assume 
that more luminous galaxies are formed by galaxy mergers 
with a shorter timescale for gas transformation into stars. 
Recently, Kobayashi (2005) and Robertson et al. (2006) have 
confirmed the importance of dissipation and the timescale 
for SF to explain mass homology breaking in ellipticals. 

In this paper we go an step further and study mass and 
velocity distributions of two samples of virtual ellipticals 
formed in self-consistent cosmological simulations. As ex- 
plained in Section 3, the simulations provide us with clues 
on the physical processes involved in elliptical formation. 
They indicate that most of the dissipation involved in the 
mass assembly of a given ELO occurs in the violent early 
phase at high z and on very short timescales (and earlier 
on and on shorter timescales as the ELO mass grows, see 
details in DSS04 and DTal06), as a consequence of ELO as- 
sembly out of gaseous material and its transformation into 
stars. Moreover, the dissipation rate history is reflected by 
the star formation rate history. During the later slower phase 
of mass assembly, ELO stellar mass growth essentially oc- 
curs without any dissipation and the SF rate substantially 
decreases (see more details and Figure 1 in DTal06) . So, the 
mass homology breaking appears in the early, violent phase 
of mass assembly and it is essentially preserved during the 
later, slower phase. A consequence is that the dynamical 
plane appears in the violent phase and is roughly preserved 
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along the slower phase, see discussion in DTal06 and Ohorbe 
et al. 2006. We see that our results on the role of dissipative 
dynamics essentially include previous ones, but they also 
add important new informations. First, our results on the 
role of dissipative dynamics to break mass homology agree 
with the previous ones, but it is important to note that, 
moreover, ELO stellar populations show age effects, that is, 
more massive objects produced in the simulations do have 
older means and narrower spreads in their stellar age distri- 
butions than less massive ones (see details DSS04); this is 
equivalent to downsizing (Cowie et al. 1996) and naturally 
appears in the simulations, so that it need not be considered 
as an additional assumption. Second, the preservation of the 
dynamical plane in the slow phase of mass aggregation in our 
simulations also agrees with previous work based on dissi- 
pationless simulations of pre-prepared mergers (Capelato et 
al. 1995; Dantas et al. 2003; Gonzalez- Garcia & van Albada 
2003; Nipoti et al. 2003; Boylan-Kochin et al. 2005). But, 
again, it is important to note that the important decrease 
of the dissipation rate in the slow phase of evolution nat- 
urally appears in the simulations and we do not have to 
assume this decrease. Also, the decrease of the merger rate 
in the later phase of mass assembly results from the global 
behaviour of the merger rate history in the particular cos- 
mological context we have considered. Third, it turns out 
that the physical processes involved in ELO formation un- 
veiled by our simulations, not only explain mass homology 
breaking (and its implications in the formation and preser- 
vation of the dynamical plane), and stellar age effects or 
downsizing in ellipticals, but they might also explain other 
elliptical properties recently inferred from observations (for 
example, the appearance of blue cores, Menanteau et al. 
2004; the increase of the stellar mass contributed by the el- 
liptical population since higher z, Bell et al. 2004; Conselice, 
Blackburne, & Papovich 2005; Faber et al. 2005; see more 
details in DTal06). 

8.4 Conclusions 

We conclude that the simulations provide an unified scenario 
where most current observations on ellipticals can be inter- 
related. In particular, this scenario proofs the importance 
of dark matter haloes in relaxed virtual ellipticals, and sug- 
gests that real ellipticals must also have extended, massive 
dark matter haloes. Also, this scenario explains the homol- 
ogy breaking in the relative dark- to bright-mass content and 
distribution of ellipticals, that could have important impli- 
cations to explain the physical origin of the Fundamental 
Plane relation and its preservation. It is worth mentioning 
that this scenario shares some characteristics of previously 
proposed scenarios, but it has also significant differences, 
mainly that most stars in EGs form out of cold gas that had 
never been shock heated at the halo virial temperature and 
then formed a disc, as the conventional recipe for galaxy 
formation propounds (see discussion in Keres et al. 2005 
and references therein). The scenario for elliptical forma- 
tion emerging from our simulations has the advantage that 
it results from simple physical laws acting on initial condi- 
tions that are realizations of power spectra consistent with 
observations of CMB anisotropies. 
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